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ABSTRACT 

Using observational data from the literature, we show that the non-thermal radio 
luminosity (L) of supernova remnants (SNRs) is a strong function of the average gas 
surface density (Sg) of the galaxy in which the remnants reside, from normal spirals to 
dense luminous starbursts. Our result supports the interpretation of the radio sources 
in M82 and Arp 220 as normal SNRs, and not "radio" supernovae. We combine a simple 
theory for electron cooling in SNRs with their observed radio luminosities to estimate 
the remnant magnetic field strength (Bsnr): the correlation between L and Eg implies 
that -BsNR also increases with Eg. We explore two interpretations of this correlation: 

(1) i?sNR is generated by post-shock magnetic field amplification, with Bg^^ oc Eg and 

(2) i?sNR results from shock-compression of the ambient interstellar medium (ISM) 
magnetic field (i?isM), with Bism being larger in denser galaxies. We find that shock 
compression is, on average, sufficient to produce the observed radio emission from 
SNRs in the densest starburst galaxies; amplification of post-shock magnetic fields is 
not required. By contrast, in normal spirals modest post-shock field amplification in 
some remnants (a factor of few — 10) is consistent with the data; we find tentative 
evidence that both the Alfven speed within SNRs and the ratio of Bsnr/^"" 
post-shock pressure ("es") are constant in SNRs from galaxy to galaxy. We discuss 
observational tests that can be used to more definitively distinguish between these two 
interpretations of the radio luminosities of SNRs. Regardless of which is correct, the 
radio emission from SNRs provides an upper limit to -Bjsm that is independent of the 
minimum energy assumption. For the densest starbursts, the magnetic energy density 
in the ISM is below the total ISM pressure required for hydrostatic equilibrium; thus 
magnetic fields are not dynamically important on the largest scales in starbursts, in 
contrast with spiral galaxies like our own. This dichotomy may have implications for 
galactic dynamo theory. 

Key words: ISM: supernova remnants — galaxies: magnetic fields, starburst — radio 
continuum: galaxies 



1 INTRODUCTION 

The magnetic energy density of the Galaxy is in rough 
equipartition with both the cosmic ray and the turbulent 
energy densities. These components of the ISM combine to 
produce a midplane pressure sufficient to balance the self- 
gravity of the Galactic disk and thus to establish hydrostatic 
equilibrium (Boulares & Cox 1990) . Moreover, direct probes 
of magnetic fields in dense star-forming regions within the 
Galaxy indicate that they may be dynamically important 



(Crutcher 1999). Because of the potential importance of 
magnetic fields for star formation on small scales, at and 
below the scale of giant molecular clouds, and for the struc- 
ture and self-regulation of galaxies on their largest scales, 
it is of considerable interest to understand how the average 
magnetic field strength varies from normal galaxies to star- 
bursts, where the ISM conditions are considerably different. 

A rough upper limit to the average magnetic field 
strength in a galaxy can be constructed by assuming that 
the field is dynamically comparable to gravity. For a thin 
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self-gravitating gas-dominated disk of average gas surface 
density Eg, this "equipartition" field strength is 

Boq ~ 27r^2GEg 2.3Eg mG, (1) 

where Eg is measured in g cm~^. Boq provides an upper limit 
to B on galactic scales because for B > Bcq the field is buoy- 
ant and escapes the host galaxy (Parker 1966). For galaxies 
like the Milky Way, in which the total surface density, Etot , 
is dominated by stars, Bcq in equation should be mul- 
tiplied by a factor of (Etot/Eg)^''^ (~ 3). The measured gas 
surface densities of galaxies vary by more than 4 decades 
(e.g., Kennicutt 1998, hereafter K98). Over this range, equa- 
tion IT]) implies a maximum field strength of Bcq ~ 15 /iG 
for the Milky Way (Eg 2 x 10"^ g cm~^; Etot/Eg ^ 10) 
and ~ 20 mG in galaxies like Arp 220 that anchor the high- 
density end of the Schmidt Law with Eg ~ lOg cm~^ and 
Etot/Eg « 1 (Downes & Solomon 1998). 

Estimates of magnetic field strengths in galaxies are 
traditionally limited to the "minimum energy" assumption 
of Burbidge (1956), which posits an equality between the 
magnetic and cosmic ray energy densities (see Longair 1994; 
Beck & Krause 2005). This approximation works well in the 
Galaxy and other normal spirals where the field strength is 
also nearly equipartition in the sense of equation ([l|, but 
the minimum energy assumption likely underestimates the 
field strengths in starburst galaxies because of strong cosmic 
ray electron cooling (Condon et al. 1991; Chi & Wolfendale 
1993; Thompson et al. 2006; hereafter [T06]). 

One way to see that the minimum energy estimate 
must fail in ultra-luminous infrared galaxies (ULIRGs) like 
Arp 220 is to note that the inverse Compton (IC) cooling 
time for radio-emitting cosmic ray electrons is tjc ~ 10* yr 
(e.g., Condon et al. 1991). The minimum energy estimate 
for the magnetic field strength yields a synchrotron cooling 
timescale at GHz frequencies that is ~ 10 times longer, and 
yet Arp 220 lies on the FIR-radio correlation together with 
essentially all star-forming galaxies (Condon 1992; Yun et 
al. 2001), which have synchrotron cooling times somewhat 
shorter than IC cooling times (e.g.. Fig. 2 of T06). These 
facts imply that the true magnetic field strength in ULIRGs 
is considerably larger than the minimum energy estimate, 
and that rapid electron cooling invalidates the assumptions 
upon which that estimate is predicated. 

This argument can be generalized and used to con- 
struct an empirically derived minimum magnetic field 
strength in galaxies by equating the magnetic energy 
density Ub and the energy density in starlight Uph- Because 
the ratio Ub /Uph determines the ratio of synchrotron to 
IC cooling, and because normal galaxies and ULIRGs lie 
on both the FIR-radio correlation and the Schmidt Law, 
Ub = B'^ /8iT must be larger than (or a constant fraction of) 
C^ph = F/c — eEsFRC, where e is a stellar IMF-dependent 
constant, F is the fiux, and Esfr is the star formation rate 
per unit areaQ Thus, the average magnetic field strength 
must exceed 



^ In principle, a factor of the dust optical depth should be applied 
to Uph in the densest starburst galaxies, if the cosmic rays are co- 
spatial with the dense molecular gas that is optically thick even in 
the FIR (see, e.g., Thompson et al. 2005). However, this correction 
is uncertain and for the purposes of constructing a lower limit to 
B based on Uph, this expression suffices. 



Bph = (87reEsFRc)^''^ ~ 0.3Eg'^mG 

^ 1.0 E;^ *'' mG, (2) 

where the numerical approximations follow from the ob- 
served Schmidt Laws of K98 (top) and Bouche et al. (2007) 
(bottom). Bph is a lower limit to the mean field strength 
in galaxies because if Ub < C^ph, variations in UB/Uph 
from galaxy to galaxy would likely introduce scatter and 
non-linearity into the FIR-radio correlation (T06). For the 
Galaxy, B > Bph ~ 4 /iG and for Arp 220, B > Bph « 1 mC. 

In this paper, we argue that radio observations of su- 
pernova remnants (SNRs) provide a complimentary probe 
of the ISM magnetic field strength in star-forming galaxies. 
In fj2] we describe a simple synchrotron cooling model for 
the radio luminosity of an individual SNR. The radio fiux 
depends critically on the magnetic field in the SNR, which 
is somewhat uncertain. At a minimum, a SNR must contain 
shock-compressed ISM magnetic field, which is constrained 
(on average) by the upper and lower limits of equations ([1} 
and ([2]), respectively. However, the magnetic field in SNRs 
could be much stronger than the shock compressed ISM field 
if there is significant amplification in the post-shock plasma. 
We compare the remnant magnetic field strengths derived 
from the model presented in fj2]with our expectations from 
equations ([T| and ((2]), to set limits on the importance of 
field amplification in SNRs. We further argue that the radio 
emission from SNRs provides an upper limit on the ambient 
ISM field in star-forming galaxies. In ^ we show that the 
luminosity of SNRs is strongly correlated with the average 
gas surface density of the galaxy in which they reside. Using 
the model of ^ we invert these observations to constrain 
the remnant and ambient ISM magnetic field in this sample 
of galaxies. We discuss our results in SJH focusing on the rel- 
ative importance of shock compression of ISM field versus 
post-shock field amplification. 



2 MAGNETIC FIELDS IN SUPERNOVA 
REMNANTS 

To estimate the radio luminosity of a SNR, we assume that 
a fraction 10"^^ of the supernova kinetic energy {E51 — 
Bsn/IO^^ ergs) is supplied to primary cosmic ray electrons 
in the shock and that the accelerated electrons radiate syn- 
chrotron in a magnetic field of strength BmG = B/mG. 
We further assume that the electron particle spectrum is 
fiat, with a power-law index of p = 2 (71(7) oc 7"*"), as is 
expected theoretically for strong shocks (e.g., Blandford & 
Eichler 1987) and is observed in situ in some SN remnants 
(e.g., Aharonian et al. 2005; Brogan et al. 2005). A value of 
^ ~ 1 is required to explain the integrated radio flux from 
star- forming galaxies, i.e., the FIR-radio correlation (Volk 
1989; T06). Given the small scatter in the FIR-radio corre- 
lation (Yun et al. 2001), f « 1 is uncertain at the factor of 
< 2 level when averaged over many SNe in a galaxy. With 
these assumptions, the radio luminosity is simply (van der 
Laan 1962) 

""^•^ ~ Tw^U, vr~ ~ 3 >^ 10 ? -^51 '^GHz -B ergs s , (3) 
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where tsyn ~ 3 x 10* bJ^^ '^ghz ^^e synchrotron cool- 

ing timescale for electrons radiating at j^ghz = i^/GHz and 
where we have assumed that tsyn > ^exp, where tcxp is the 
expansion time of the remnant. In the last estimate in equa- 
tion ((3]) we have assumed 7max = 10®; variations in 7max 
from 10* — 10* affect this estimate by only a factor of ~ 1.5. 

Equation Q is an estimate of the radio emission pro- 
duced by the interaction between a SN and the ambient 
ISM, i.e., the emission from a SNR (e.g., Shlovskii 1960). 
The radio emission is expected to peak at the Sedov time 
and persist until the remnant cools radiatively (see the dis- 
cussion of these timescales in JH). This time evolution is in 
contrast to radio supernovae which peak and decay after a 
month to few-years and which likely result from the inter- 
action between a SN shock and the immediate circumstellar 
medium and pre-SN ejecta (e.g., Chevalier 1982; Weiler et 
al. 2002). Nearly all of the radio remnants discovered in M82 
and Arp 220 that we discuss below do not vary apprecia- 
bly on decade timescales (Kronberg et al. 2000; Rovilos et 
al. 2005) . Rather than interpret these bright sources as bona 
fide radio SNe, we thus interpret them as normal SNRs. 

Given an observed radio flux density, equation (|3} can 
be inverted to estimate the magnetic field strength in the 
SNR, -BsNR. Assuming that Bsnr is ambient ISM field com- 
pressed by the SN shock by a factor of / « 3 (appropriate 
for randomly-oriented fields; Volk et al. 2002 advocate / ~ 6 
for young remnants), we obtain a simple estimate of the ISM 
magnetic field strength in terms of the SNR flux density: 

equivalent to Bism/ ~ Bsnr- Equation Q is most readily 
interpreted as an upper limit on the ISM magnetic field given 
a set of SNRs with flux density Si, because (1) / may be 
as large as ~ 7 and (2) turbulence and instabilities can in 
principle amplify the post-shock magnetic field. The SNR 
field may thus be larger than ~ /-Bism , but it is unlikely to 
be smallerQ 

The assumptions made in deriving equation (|4]) include 
that tsyn > texp (wcak cooling) and that p — 2. Together, 
these assumptions imply that the synchrotron spectral index 
should be a = (p l)/2 = 1/2 (S^ oc i/""). For p / 2, 
and fixed total energy injected in cosmic ray electrons, the 
normalization of equation Q increases, but modestly: for 
p = 2.6, -Bism is larger by a factor of « 1.5 and iJisM oc S^^^ . 
The spectral indices of GHz radio emission from SNRs range 
from a ~ 0.5—1 (see the refs. in Table [l|. Diversity in 
a may signal diversity in p, evolutionary effects, that the 
assumption of weak cooling is violated, or some combination. 
Additionally, note also that ^ ~ 1 is not valid during the 
early free-expansion phase of SNR evolution. 

A further assumption made in deriving equation (|4]) is 
that the supernova remnants are not radiative. If they were, 
/ could be arbitrarily high. Indeed, Chevalier & Fransson 
(2001) argued that the radio point sources observed in M82 
and Arp 220 are radiative, but otherwise normal, supernova 
remnants (i.e., not "radio" SNe). As we show below and as 

If SNe systematically sample only low magnetic field regions of 
a galaxy, -Bism inferred from SNRs may not be the same as the 
volume averaged magnetic field strength. This possibility cannot 
be ruled out, but we regard it as unlikely. 




log,o[t (yr)] 

Figure 1. The time evolution of the post-shock magnetic field 
strength, assuming -B^/Stt = esPshi where Pgi^ is the post-shock 
thermal pressure. A 10^^ ergs supernova shock and £b = 0.1 
are assumed in all cases. With these assumptions, the post-shock 
magnetic field is strongest during the free expansion phase, and 
declines during the Sedov- Taylor and snow-plow phases. 



we discuss in more detail in 21 the radio luminosities of 
the remnants in starburst galaxies are consistent with com- 
pression of the ambient ISM magnetic field by an adiabatic 
shock, perhaps together with modest post-shock fleld ampli- 
flcation. Much larger /, as would be implied if the remnants 
were radiative, is inconsistent with independent constraints 
on the ambient ISM magnetic fields in these systems. 

To provide some quantitative context for comparing the 
strength of the magnetic field in SNRs produced by fiux 
freezing with that produced by possible amplification in the 
post-shock plasma. Figure [T] shows the post-shock magnetic 
fleld strength in SNRs as a function of time assuming that 
the post-shock remnant magnetic energy density is a fac- 
tor es times the post-shock pressure (as is often assumed in 
models of radio supernovae and gamma-ray burst afterglows; 
e.g., Reynolds & Chevalier 1981); we consider a range of am- 
bient densities n, from low densities appropriate to SNRs 
in the Milky Way to high densities appropriate to SNRs in 
Arp 220. For these calculations we have used the analytic ap- 
proximations for the evolution of SNRs developed by Draine 
& Woods (1991), which include both the Sedov- Taylor and 
pressure-driven snow-plow phases. The post-shock magnetic 
fleld is the strongest during the free expansion phase and 
decreases during the Sedov- Taylor phase and later as the 
post-shock plasma pressure decreases. We note, however, 
that the appropriate value for es is quite uncertain, and 
thus it is difficult to estimate from flrst principles the im- 
portance of post-shock field amplification (see Riquelme & 
Spitkovsky 2009 for a recent study). 
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Figure 2. SNR radio luminosity versus average gas surface density, Sg, for the sample of star-forming galaxies in Table [T] The filled 
square and horizontal bar show the maximum and average SNR luminosity for each system. The vertical bar extends to the lowest SNR 
luminosity observed. A strong correlation is present in the data. A simple linear least squares fit to the average SNR luminosities gives 
uLi, oc S^-*^; fitting to only the brightest remnants gives i/L^ oc (see eq.[5]l. Note that although for some galaxies many SNRs 

are observed (e.g., M33 with 5f SNRs), others, like 2146, have only a handful, or in the case of Arp 299, a single confirmed SNR. For a 
similar figure, see Hunt & Reynolds (2006). 



3 RESULTS 

A large sample of radio SNRs is available in the literature. 
These SNRs reside in a variety of galaxies that span the full 
range of both the FIR-radio correlation and the Schmidt 
Law. In general, the remnants are identified using the VLA, 
MERLIN, or VLBI, and are distinguished from compact HII 
regions by their steep non-thermal spectra {a > 0.5). The 
few "radio SNe" proper that appear in the compilations of 
SNe presented in the references in Table [T] are excluded. 
These include remnant 41.954-57.5 in M82 from Kronberg et 
al. (1985) (see Kronberg et al. 2000) and SN 2000ft in NGC 
7469 (Colina et al. 2001; Alberdi et al. 2006). Background 
objects such as AGN behind nearby galaxies like NGC 5194 
are a contaminant to the SNR sample (see Table [!}. 

Figure [2] shows the luminosity of observed SNRs as a 



function of the average gas surface density Eg of the galaxy 
in which they reside. For simplicity, all fluxes have been 
scaled to a frequency of 1.4 GHz assuming a = 0.5, even in 
those cases where a is known empirically in some frequency 
range. References, assumed distance, number of remnants, 
maximum and average SNR luminosity, and the sensitiv- 
ity, frequency, and resolution of the observations are given 
in Table [T] For all of the galaxies we consider, the remnants 
account for just a few to ten percent of the total non-thermal 
continuum. Each vertical line in Figure[2]is a separate galaxy 
(or, for Arp 220, component). The filled square shows i/L,j 
for the brightest single SNR. The line for each system ex- 
tends to the lowest luminosity remnant identified. The hor- 
izontal line shows the average luminosity of the observed 
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Figure 3. Inferred ISM magnetic field strength as function of Eg for the galaxies in Table [T] (from eq.|4]with / = 3), assuming that shock 
compression of ambient ISM field determines the magnetic field strength in SNRs (Ssnr); alternatively, this plot can be interpreted as 
showing -BsNR.//- As in Figurc[2] the filled square and horizontal bar show the maximum and average value of B from the brightest and 
average SNR, for each system. The vertical bar extends to the lowest SNR luminosity observed. The thick solid line is the maximum 
magnetic field estimate for a gas-dominated self-gravitating disk, Bcq (eq. [I]). The dotted line includes the correction discussed after 
eq. JTJl for disks that are not gas-dominated (Stot/Sg = 10). The dashed lines show the minimum required magnetic field strength for 
consistency with the FIR-radio correlation (see eq. [2]l, assuming the average Schmidt Law from K98 (shallower index) and Bouche et 
al. (2007) (steeper index). An unweighted linear least-squares fit to -Bism or ^SNRi inferred from the average luminosity SNR (horizontal 
lines on each bar), gives B oc Sg'^^, whereas a fit to the brightest SNRs (highest inferred field strengths, filled squares) gives B oc Sg'^^ 
(compare with Fig. |4jl. 



SNR distribution. A similar figure was constructed by Hunt 
& Reynolds (2006). 

A strong correlation is evident between i/Li, and Eg. A 
simple power-law fit to the average SNR luminosities (hori- 
zontal bars) gives 



logio 



^L^(1.4GHz) 
ergs s^i 



36.03 + 0.85 log 



g cm^ 



(5) 



A similar fit is found for the brightest SNRs (filled squares), 
but with a normalization of 36.48 and a slope of 0.78. It is 
also clear that the lowest SNR luminosities exhibit a cor- 



relation with Eg that is qualitatively similar to that of the 
average and brightest SNRs. 

The latter raises the concern that the results of Figure 
[2] may represent a selection bias rather than a physical ef- 
fect. The correlation between the maximum SNR luminosity 
and Eg in Figure [2] is the most secure and is not subject to 
an obvious bias. In particular, although faint remnants in 
high-Eg galaxies (which are rarer and thus tend to be more 
distant) would be undetectable, bright remnants could be 
readily seen in nearby low-Eg systems, but are not. On the 
other hand, the correlations between the average and lowest 
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remnant luminosity and Eg are more subject to the obvious 
bias that faint remnants are difRcult to detect in all galaxies. 
One might anticipate that this bias would be particularly se- 
vere for higher Eg galaxies and could account for the absence 
of low luminosity remnants in these systems. One indication 
that this effect does not significantly affect our calculation 
of the average SNR luminosity is that nearly all of the de- 
tected remnants in M82, at a distance of 3.6 Mpc, are more 
luminous than even the brightest remnant in NGC 300 at 
a distance of 1.9 Mpc (see Table [!}. An additional piece of 
information comes from comparing the samples of Muxlow 
et al. (1994) with Fenech et al. (2008) for the remnants in 
M82. The latter are three times more sensitive than the for- 
mer, and yet, comparing the two datasets we see that the 
average SNR luminosity decreases by just a factor of ~ 1.5, 
from log 10 ['^i™'"'"] ~ 35.1 to fa 34.9. A detailed comparison 
of the luminosity function of detected SNRs in all galaxies 
in this sample would be useful both for exploring the den- 
sity and magnetic field distribution of the ISM ( i)4.3|l . and 
for a more detailed understanding of the selection effects at 
low remnant luminosity. Finally, we note that an additional 
reason that the lowest luminosity SNRs correlate with Eg is 
that detectable remnants must be brighter than the diffuse 
radio continuum of their host galaxy, which is larger in high- 
Eg galaxies as a result of the Schmidt Law (Esfr oc Eg *; 
K98; Esfr oc EJ ''; Bouche et al. 2007) and the FIR-radio 
correlation. 

From the observed luminosities in Figure (2] we can in- 
fer the SNR magnetic field Bsnr under the assumptions of 
equations ([Sjl and Q. To compare this magnetic field es- 
timate with other constraints on the ISM magnetic field, 
we plot these results in Figure [3] in terms of the ISM mag- 
netic field strength Bism = Bsnr//, i.e., the ambient ISM 
field strength that would produce a SNR field strength of 
BsNR by flux freezing. The correlation in Figure [2] combined 
with equation Q implies that Bsnr and Bism are strong 
functions of Eg . A simple linear least-squares fit to the data 
points for the average SNRs (horizontal bars) gives 



logio 



-BiSM -Bsnr/3 



-3.20 0.565 log 1 



g cm^ 



.(6) 



For the brightest SNRs the relation is similar, but with 
normalization —2.90 and slope 0.52. Note that for larger 
assumed p, the correlation between Bism and Eg flattens 
(see the discussion after eq. U} . Figure |3] demonstrates that 
for normal spirals the magnetic fleld strength Bism inferred 
from the average-brightness SNR is comparable to or larger 
than the equipartition fleld, B^q (dotted line), and larger 
than the minimum magnetic field strength Bph (dashed line; 
eq. [2}; this is true up to Eg ~ 0.1 g cm~^, which character- 
izes starbursts such as NGC 253 and M82. For the densest, 
most luminous starbursts in our sample (Arp 299 and the 
nuclei of Arp 220) , Bjsm is below Bcq and — at least for Arp 
220 — comparable to Bph. At the other extreme, for the low- 
est density systems in our sample, NGC 300 and NGC 4736, 
Figure |3] tentatively implies a trend of larger BisM/Bcq at 

low Eg. 

Because of the paucity of systems at very high Eg in 
Figures [2] and O we searched the literature for additional 
high resolution radio observations of local starbursts and 
ULIRGs. Unfortunately, the existing observations are not 
yet at sufficiently high spatial resolution to unambiguously 
imply that single SNRs have been detected. For example. 



Momjian et al. (2006) found 10 compact radio sources in 
the ULIRG IRAS 17208-0014. Most have deconvolved sizes 
at half maximum of ~ 10 pc. At a distance of ~ 170 Mpc 
{z — 0.043), it is unclear if these compact sources are 
indeed individual SNRs, rather than unresolved popula- 
tions of multiple remnants. In addition, the spectral indices 
of these sources are unknown, and many may be bright 
HII regions. A single source is unresolved, with a linear 
scale < 4.9 pc. This is also the least luminous source with 
5^(1.6 GHz) — 119 ± 45 ^Jy. Given the gas surface density 
from K98 (~ 0.4 g cm~^), this individual source would lie 
above the single Arp 299 remnant (see Table [l| in Figure [2] 
by a factor of roughly two. Similarly, recent observations of 
IRAS 23365-H3604 and IRAS 07251-0248 have been reported 
by Romero-Canizales et al. (2008). They reveal a number of 
compact radio components. However, because these galaxies 
are at distances of « 230 and « 360 Mpc, respectively, the 
resolution of the observations is again insufficient to asso- 
ciate individual sources with single SNRs. Comparing with 
the observations of the nuclear sources in Arp 220 by Lons- 
dale et al. (2006), we expect that the sources detected by 
Romero-Cafiizales et al. (2008) are likely to be collections of 
~ 10 SNRs per source. Higher resolution studies can confirm 
or exclude this possibility. Taking the sources from Romero- 
Canizales et al. (2008) at face value, and estimating the gas 
surface density of these ULIRGs from the literature, both 
systems would lie significantly above the solid line (Beq) in 
Figure |3l and they would deviate markedly from the trend 
seen in Figure [51 



4 DISCUSSION 

Figure [2] demonstrates that SNR luminosities are a strong, 
and fairly continuous, function of Eg, the average (galaxy- 
wide) gas surface density of their host galaxy. There may be 
a change in the correlation for Eg > 0.1 g cm~^, but more ob- 
servations of SNRs in dense starbursts are required to assess 
this. In the following, we consider two plausible interpreta- 
tions of this data: (1) the luminosity of a SNR is significantly 
larger in denser galaxies because the ambient magnetic field 
is significantly larger (as in Fig. [3)l , or (2) SNRs are brighter 
in denser galaxies because the magnetic field in the post- 
shock plasma scales as /Siv cc pisM^sh is thus larger 
in galaxies with a denser ISM (larger pism). 



4.1 Field Amplification 

Recall that Bph (eq.[2l dashed lines in Fig. [3]) is a lower limit 
to the true average ISM magnetic field strength because for 
B < Bph IC losses dominate synchrotron losses and one ex- 
pects a systematic deviation from the FIR-radio correlation 
at high FIR luminosities, which is not observed (T06). In ad- 
dition, the magnetic field inferred from SNRs, Bism (eq. HI 
points in Fig. [3]) is an approximate upper limit to B because 
post-shock turbulence and instabilities can amplify magnetic 
fields in SNRs (Volk et al. 2002; Berezhko et al. 2006). It 
is striking that the upper and lower limits on Bism coincide 
to within a factor of few for high Eg starburst galaxies, par- 
ticularly given the maximum equipartition magnetic fields 
these galaxies could support (Boq, eq. [T]). This correspon- 
dence implies that shock compression of Bism is sufficient 
to power SNR radio emission in dense starburst galaxies. 
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(g cm-=) (g cm-^) 

Figure 4. Left Panel: SNR Alfven speed (Va) as a function of tlie average ISM gas surface density Eg. Here, pisM = ^g/C^h), where 
h = 100 pc is the assumed gas scale height for all galaxies. Note that the range of Eg plotted here corresponds to the range of average 
densities 0.5 cm~^ < njgM < 10'* cm~'^. Right Panel: The ratio of the post-shock magnetic field energy density to the post-shock thermal 
pressure, eg = [Bgj^jj^/(87r)]/[(3/4)pigM V^'^], as a function of Eg for h = 100 pc. In both panels, as in Figure|3] BsNR is inferred from 
equation Q and the SNR radio luminosities in Table [Tl Note that if supernovae explode in regions of lower-than-avorage density ( ij4.1l 
i|4.3l l. the inferred Va and eg both increase. 



For the specific case of Arp 220, our inferred ISM mag- 
netic field strength of ~ 1 — 3 mG is also reasonably con- 
sistent with the direct Zeeman detections of Robishaw et 
al. (2008); thus three independent observational methods, 
each likely probing somewhat different phases of the ISM, 
all imply that a magnetic field of a few mG pervades the 
nuclear region of Arp 220. A corollary of this is that Figure 
[3]provides strong evidence for the interpretation of the radio 
point sources in Arp 220 as normal SNRs expanding into a 
highly-magnetized ISM; this alleviates the need for bona fide 
radio SNe to generate the very high luminosities observed 
(Smith et al. 1998; Rovilos et al. 2005; Parra et al. 2007; see 
also Chevalier & Fransson 2001). 

In contrast to starbursts like Arp 220, in normal spirals 
independent observations constrain the ambient ISM mag- 
netic field to be ~ 5 - 10 ^G (e.g.. Beck 1982; Fitt & Alexan- 
der 1993; Beck & Krause 2005). Compressing this field by 
a factor of ~ 3 — 7 results in a magnetic field of sufficient 
strength to explain the typical radio SNR observed in nor- 
mal spirals (Fig. [3]); it appears, however, that the brightest 
remnants require fields that are a factor of few larger (and 
perhaps more at the lowest surface densities), either because 
the SNe explode in a region of larger-than-average B or be- 
cause of modest post-shock field amplification. Nonetheless, 
as in the case of starbursts. Figure [3] implies that significant 
post-shock magnetic field amplification (say > x 10) is not 
required to explain the average radio luminosities of SNRs 
in spiral galaxies. 

In our own Galaxy, there is some evidence for magnetic 
field amplification in young shell- type SNRs (e.g., Volk et 
al. 2002; Volk et al. 2005; Berezhlco et al. 2006). In these 



models, individual Galactic SNRs are observed and modeled 
at a moment in their evolution, in contrast to the statisti- 
cal approach taken here. For example, Volk et al. (2005) 
constrain -Bsnr ~ 300 /xG in Tycho's SNR by fitting the 
broadband X-ray and radio emission simultaneously. Using 
^ = 1 in equation Q , and using only the radio data, we find 
that B in the remnant is ~ 9/iGl3 We have also compared 
the results from Volk et al. (2005) with our calculation for 
the SNRs RCW 86, SN 1006, and Cass A. We find that 
Bsnr ~ 3-Bism ~ 9, 4, and 180fj,G, whereas Volk et al. 
find « 100, 160, 500 /iG, respectively. The difference in the 
inferred SNR field between Volk et al. and this work lies 
partially in the assumptions about the energy in relativistic 
electrons, encapsulated in our parameter ^ in equations ^ 
and Q. In order to not over-produce the radio luminosity 
with such a large B, Volk et al. (2006) decrease the cosmic- 
ray electron-to-proton ratio, effectively using a small value 
of f . In contrast, the normalization of the FIR-radio corre- 
lation, which we use here to normalize the emission from all 
individual SNRs, requires that ^ ~ 1 when averaged over a 
suitably large number of remnants. We expect that shock 
accelerated electrons will have a total energy equivalent to 
^ ~ 1 by the onset of the Sedov- Taylor phase, when the 
shocked, swept-up mass is comparable to the ejecta mass. If 
shock acceleration only becomes efficient very late in a rem- 
nant's evolution — or if the electrons rapidly diffuse out of 
the remnant on a timescale comparable to the Sedov- Taylor 
time — then the assumption of ^ ~ 1 may be unreasonable 

^ This is for an adopted distance to the Tycho SNR, of 2.4 kpc 
(from the Green catalog; see Table [TJ . 



© RAS. MNRAS OOO.ITlfTTI 



8 Thompson, Quataert, & Murray 



for a typical remnant. Although more definitive calculations 
are required to assess these possibilities, we regard both of 
them as unlikely and thus conclude that the majority of ob- 
served remnants can be interpreted assuming ^ ~ 1. In this 
case, although individual remnants, at particular moments 
in time, may have fields larger than that due to compression 
alone. Figure [3] suggests that this is not a significant effect 
for SNRs in starbursts, or for the average SNR in most nor- 
mal spirals. 

A SNR could in principle generate magnetic fields in ap- 
proximate equipartition with the post-shock thermal pres- 
sure of the gas: Bg-^j^/Sn ~ pisM^sh where pisM is the den- 
sity of the ISM in which the SN explodes and v^h is the veloc- 
ity of the SN shock at a given time (this assumption is often 
used in radio SNe and 7-ray burst models; e.g.. Chevalier 
1982; Meszaros & Rees 1997); see Figure [1] for estimates of 
the magnetic field strengths in SNRs if significant post-shock 
field amplification indeed occurs. In this interpretation, the 

higher remnant luminosities in denser galaxies (Fig. [2| are 

1/2 

due to the post-shock scaling Bsnr oc Pjg^^, rather than 

anything about the magnetic field in the ambient ISM. 

Adopting the interpretation that post-shock magnetic 

field amplification is solely responsible for generating Bsnr 

and the observed radio emission, Figures[2]and[3]have several 

interesting implications. Taking a constant gas scale height 

of ft ~ 100 pc for all galaxies, so that the average ISM density 

is pisM = nisMrrip = Eg/(2/i), Figure |3] implies that Bsnr oc 
1/2 

PiSM' there is a constant characteristic Alfven speed in 
all observed remnantsQ The left panel of Figure|4]shows this 
explicitly. Here, we have plotted Va = SsNR/(47r4pisM)^'^'^ 
as a function of Eg, assuming a constant gas scale height. A 
fit to the average SNRs gives 



logi 



Va 



km s" 



1.87 + 0.065 logi 



g cm- 



(7) 



For the brightest remnants, the normalization is higher 
(« 2.17) and the slope is flatter 0.021). In Figure H Va 
is calculated assuming that the SNR density is 4pisM, ap- 
propriate for a strong adiabatic (7 = 5/3) shock expanding 
into the ISM. If, as we discuss in ^4.31 supernovae tend to 
sample low density regions of the ISM, the implied Va would 
increase, probably by a factor of ~ 3 — 10. Despite uncer- 
tainty in the overall normalization for Va, the left panel of 
Figure |4] shows that the Alfven speed within remnants is 
roughly the same for all observed SNRs in all galaxies; if 
magnetic field amplification is at work in SNRs, this finding 
is an important clue to the physics. 

This result can equivalently be presented as a constraint 
on es, the ratio of the post-shock magnetic energy density 
to the post-shock pressure (see Fig. [ij. The right panel of 
Figure |4] shows our determination of es for the SNRs in our 
sample as a function of the average ISM gas density. For the 
average SNRs, we find that 



logio [es] « -3.04 + 0.131 logio 



g ern- 



es) 



For the brightest remnants, the normalization is consider- 

^ Although the assumption of constant gas scale height across 
all galaxies in the sample is a reasonable first approximation, they 
may have systematically smaller h at higher Eg, which would tend 
to flatten the function VA{n-isM) (see eq. (Tjl. 



ably higher —2.44) and the slope is yet flatter (« 0.042). 
The value of es is proportional to V^, so the near-constancy 
of e_B with TiisM is not surprising given the results for Va in 
the left panel of Figured As with Va, we emphasize that 
the precise value of es depends sensitively on whether or 
not supernovae explode in gas at the average ISM density. 
If, as is more likely, most of the volume is filled with gas 10 
to 100 times less dense than nisM ( ^4.3p . the typical inferred 
value of es would be ~ 0.01 — 0.1. 

Because the observations do not definitively favor either 
SN shock compression of ISM field or post-shock field ampli- 
fication, we are unable to unambiguously determine which 
of these two possibilities is correct. An important testable 
difference between these two hypotheses lies in the depen- 
dence of the radio emission on shock velocity or SNR di- 
ameter -DsNR- For shock compression alone, i/L,^ should be 
independent of Vah while if B^/Stv cx pisMWshi ^Li, oc 
Another way to state this result is that the two models 
predict different SNR surface brightness-to-diameter rela- 
tions (the EsNR — DsNK relation). For shock compression 
alone, Esnr oc -DgNR while for B oc v^^, Esnr oc -DgNR 
(in both cases assuming p = 2; see also, e.g., Berezhko & 
Volk 2004). Observations of SNRs in the Milky Way find 
Esnr oc -Dsnr (Case & Bhattacharya 1998), much closer to 
the fiux-freezing prediction, while observations of SNRs in 
M82 find SsNR oc D"^ '^ (Huang et al. 1994), closer to the 
predictions of the field amplification model (see Urosevic et 
al. 2005). For M82, the recent work by Fenech et al. (2008) 
finds a best fit of Esnr oc D""^". Updated observational 
samples, particularly those that characterize SNRs in high- 
density regions, and more sophisticated theoretical mod- 
els are needed to fully understand the implications of the 
Esnr — -Dsnr relation and its connection to the correlation 
presented in Figure [2] 



4.2 The Dynamical Role of Magnetic Fields 

Another important conclusion from Figure|3]is that Bism ^ 
Seq in the most luminous starbursts. This implies that mag- 
netic fields are not dynamically important for hydrostatic 
equilibrium on large scales in these systems, in sharp con- 
trast with spiral galaxies like our own. One possible expla- 
nation for this difference is that the lifetime of intense star- 
bursts (~ 10'' - 10** yrs) is significantly shorter than that 
of the relatively continuous star formation in normal spiral 
galaxies; there may thus be less time for dynamo processes 
to amplify the average ISM magnetic field. 

Figure [3] also suggests that the magnetic energy den- 
sity Ub is comparable to the photon energy density Uph in 
a wide range of systems, with Ub ~ few— lOf/ph, depend- 
ing on the form of the Schmidt Law adopted (compare the 
two dashed lines in Fig. [3]). This result is not as surpris- 
ing as it might first appear because massive stars determine 
both the radiation field of galaxies and the energy and mo- 
mentum injected into their ISM. The latter "feedback" can 
generate turbulence in the ISM, amplifying the magnetic 
field. Thompson et al. (2005) and Thompson (2008) have 
suggested a direct connection between turbulence in dense 
starbursts and Uph by positing that radiation pressure from 
the absorption and scattering of stellar light by dust grains 
dominates the overall pressure support. The rapid diffusion 
of infrared photons in starbursts causes the medium to be 
gravitationally unstable (Thompson 2008) and, even in the 
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absence of self-gravity, generates magnetohydrodynamic in- 
stabilities that can amplify magnetic fields directly (Turner 
et al. 2007). In addition, the momentum supplied to the ISM 
by stellar winds and SN explosions is comparable to that in 
photons. There are thus a number of sources capable of gen- 
erating turbulent motions that could amplify magnetic fields 
to the levels inferred here. It is also worth noting that al- 
though Bph <C -Beq for luminous starbursts in Figure O this 
does not necessarily imply that radiation pressure is dynam- 
ically unimportant; the dense gas in starbursts is optically 
thick even in the far infrared and thus the energy density of 
photons in the densest gas is larger than that in equation 
((21) by a factor of the dust optical depth (see footnote I), 
which can be ~ 10 — 100 (e.g., Thompson et al. 2005). 

4.3 SNR Luminosity Function 

As a final comment, we note that a more sophisticated un- 
derstanding of the SNR luminosity function is needed to 
provide a more quantitative interpretation of the data in Fig- 
ures[2]&|3] For example, although equation ([3]) suggests that 
vLi, is independent of time, the radio luminosity of a given 
SNR will peak at its Sedov time and then decline when the 
remnant cools and the shock loses energy (the radio luminos- 
ity can decline more quickly if tsyn is short or if the relativis- 
tic electrons rapidly diffuse out of the remnant; see, e.g., van 
der Laan 1962; Baring et al. 1999; Berezhko & Volk 2004). If 
the average ISM magnetic field is volume-filling, and if SNe 
sample the medium fairly, taking an average density of Arp 
220 of (n) ~ 10'' cm"^ and a SN rate of ~ lyr"^ we esti- 
mate that ~ 10^ SNRs should be visible in the radio, which 
is at odds with the ~ 50 SNRs identified (e.g., Lonsdale et 
al. 2006). However, there is an important set of mitigating 
effects. First, calculations of high Mach number isothermal 
turbulence imply that for the conditions of Arp 220 the aver- 
age medium is filled with gas with n ~ 10^^ (n) ~ lO'^ cm~'^ 
(e.g., LeMaster & Stone 2008). Second, although the syn- 
chrotron cooling time tsyn is longer for lower ambient den- 
sity (because of an assumed smaller -Bism), and one might 
thus expect to see more remnants, the SNRs at t ~ tsyn are 
also physically larger: for n — 10^ cm~^ a remnant expands 
to 4pc by ~ tsyn. For Arp 220, this is large enough that 
the SNRs would overlap on the sky and be resolved out by 
current VLBI observations. The radio luminosity function 
of SNRs is thus a convolution of the spatial variation of the 
ISM conditions into which SNe explode (e.g., variations in 
PiSM and -Bism), the time- dependence of an individual rem- 
nant, and the sensitivity and resolution of the observation. A 
careful study of these effects might provide important con- 
straints on the nature of the ISM — the statistics of the 
density field and its magnetization — in external galaxies. 
In addition, it would more quantitatively determine the im- 
portance of post-shock magnetic field amplification in SNRs, 
and how the average (or maximum) SNR luminosity is re- 
lated to the host galaxy's volume- averaged ISM magnetic 
field strength. 
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Table 1. Radio Emission from Supernova Remnants 



System A^snr'' ^g" log[!/L™<=='"] logfi/L™^"] = Freq.-f BeamS Sensitivity'' Refs. Refs. 

Name (Mpc) (gcm-^) (1.4 GHz) (1.4 GHz) (GHz) SNRs Eg 

(ergs s~-^) (ergs s"-') 



NGC 300 


i.gp 


17 


0.00050 


33.4 


33.7 


1.45 


4.7" X 3.6" 


eO/xJy 


1 


2 


NGC 4736 


4.4 


10 


0.00093 


34.3 


34.5 


1.45 


... 5 


... "J 


3 


3 


LMC 


0.053 


21 


0.0017' 


33.6 


34.4 


0.48 


4'. 3 


few m Jy 


4,5 


6 


Milky Way 


n 


38° 


0.0020 


33.2 


34.5 








7 


8 


M33 


0.84 


51 


0.0022 


33.1 


33.8 


1.42 


7" 


50AtJy 


9 


10 


SMC 


0.065 


23 


0.0033^ 


32.9 


33.5 


2.37 


40" 


0.4 mjy 


11 


6 


NGC 6946 


5.6 


15 


0.0042 


34.3 


34.6 


1.4 


2" 


20MJy 


12 


10 


NGC 1569 


2.2" 


23 


0.0044 


33.8 


34.6 


1.49 


1.4" X 1.4" 


21/xJy 


13 


10 


M51 


6.5 


37 


0.0062 


34.4 


35.2 


1.4 


1.50" X 1.21" 


22.5 AtJy 


14 


10 


M83 


7.3 


4 


0.011 


34.7 


34.9 


1.4 


3.5" X 3.5"'' 


45/iJy'' 


15 


10 


NGC 253 


3.2 


11 


0.067'^ 


35.3 


35.9 


5.0 


0.6" X 0.3" 


40MJy 


16 


10,17 


NGC 2146 


14.5 


3 


0.076 


36.1 


36.3 


1.6 


0.19" X 0.15" 


35 fj.]j 


18 


10,19 


M82 


3.6 


37 


0.081' 


34.9 


35.9 


5.0 


35 — 50 mas 


17AtJy 


20 


21 


NGC 4945 


3.8 


8 


0.19™- 


35.6 


36.0 


2.3 


16 — 15 mas 


75 AtJy 


22 


23 


Arp 299 


41 


1 


2.6 


36.6 


36.6 




5.6 mas X 4.5 mas 


35 fj.]j 


24 


10 


Arp 220E 


78 


20 


6.0 


36.1 


36.5 


1.65 


5.9 mas X 2.7 mas 


5 — 6 fj.]y 


25 


26 


Arp 220W 


78 


29 


8.7 


36.6 


37.1 


1.65 


5.9 mas X 2.7 mas 


7-9AtJy 


25 


26 



"Distance in Mpc. 

''Number of supernova remnants. 

■^Average gas surface density. 

''Mean SNR luminosity. Each SNR luminosity is calculated assuming oc u~'^'^. 
"^Maximum SNR luminosity. 

^Frequency of quoted flux density in SNR reference used to calculate SNR, luminosity. 
^Measure of beamsize quoted in SNR reference. 

''Measure of sensitivity quoted in SNR reference. Typically RMS sensitivity per beam. 
'Calculated using a total gas mass of 6 X 10* Mq (Israel 1997) and R25 ~ 4.9 kpc. 
^Calculated using a total gas mass of 4.5 X 10* Mq (Israel 1997) and R25 ~ 3.0 kpc. 

''Surface density scaled to K98, but with the Xqq conversion factor advocated in Mauersberger et al. (1996a). 
'Total gas mass of starburst region 2.3 X 10* Mq (Weiss et al. 2001), with diameter ^ 870 pc, adjusted for D = 3.6 Mpc. 
""Total gas mass within a radius of 12"(ai 221 pc at D = 3.8 Mpc) is taken as 1.5 X 10* Mq (Mauersberger et al. 1996b). 
"All distances from Green (2006); see also Case & Bhattacharya (1998). 

°Prom Green catalog. Only "S"-type sources with unambiguous distances are included. For comparison, using a sub-sample of the 
20 SNRs from Case & Bhattacharya (1998) that overlap with the Green catalog, we have compared the highest and average SNR 
luminosities using distances from the former and the latter, respectively. Both the highest SNR luminosity (Cass A) and the average 
SNR luminosity are unchanged: logjQ [i'L™''™(1.4 GHz; ergs s~^)] ^ 34.5 and ~ 33.5, respectively. However, for the full Green catalog of 
38 sources (listed above), logj^Q[!^L™'^^"(1.4 GHz; ergs s~^)] decreases to 33.2, a factor of Kb 2.0 lower than the smaller sample of Case & 
Bhattacharya (1998). 
^Distance from Gieren et al. (2005). 
''Values not listed in reference. 

""See Maddox et al. (2006), their Table 1. The four sources listed are those with confirmed optical counterparts. 
"Distance from Israel (1988). 

References: (1) Pannuti et al. (2000); see also Payne et al. (2004); (2) Read et al. (1997); (3) Duric & Dittmar (1988); 
(4) Mathewson et al. (1983); (5) Mathewson et al. (1984); (6) Israel (1997); (7) fluxes from Green, D. (2006) (see 
http://www.nirao.cam.ac.uk/surv6ys/snrs/); see also Case & Bhattacharya (1998); (8) Boulares & Cox (1990); (9) Gordon et al. (1999); 
(10) K98; (11) Filipovic ot al. (2005); (12) Hyman et al. (2000); (13) Chomiuk & Wilcots (2009); see also Greve et al. (2002); (14) Maddox 
et al. (2007), including all sources steeper than a = 0.4, excluding sources 53 (the nucleus) and 104 (background; Ho & Ulvestad 2001); 
(15) Maddox et al. (2006); (16) Ulvestad & Antonucci (1997), see also Ulvestad & Antonucci (1991); (17) Mauersberger et al. (1996a); 
(18) Tarchi et al. (2000); (19) Greve et al. (2006); (20) Fenech et al. (2008), but see also Kronberg et al. (1985), Huang et al. (1994), 
Muxlow et al. (1994); Wills et al. (1997), and Allen & Kronberg (1998); (21) Weiss et al. (2001); (22) Lenc & Tingay (2008) (23) 
Mauersberger et al. (1996b); (24) Neff et al. (2004); (25) Lonsdale et al. (2006), and see also Rovilos et al. (2005), Parra et al. (2007); 
(26) Downes & Solomon (1998) 
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